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In supersymmetric models a gluino can decay into tbxi through a stop or a sbottom. The decay 
chain produces an edge structure in the m t b distribution. Monte Carlo simulation studies show 
that the end point and the edge height would be measured at the CERN LHC by using a sideband 
subtraction technique. The stop and sbottom masses as well as their decay branching ratios are 
constrained by the measurement. We study interpretations of the measurement in the minimal 
supergravity model. We also study the gluino decay into tbxf as well as the influence of the stop 
left-right mixing on the nibb distribution of the tagged tb events. 
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I. INTRODUCTION 

The minimal supersymmetric standard model (MSSM) 
is one of the promising extensions of the standard model 
[1]. The model requires superpartners of the standard 
model particles (sparticles), and the large hadron collider 
(LHC) at CERN might confirm the existence of the new 
particles [2]. The LHC is a pp collider at ^/s = 14 TcV, 
whose operation is currently expected to start in 2007. 
The integrated luminosity will be 10 fb _1 /year at the 
beginning (low luminosity runs), and then upgraded to 
100 fb /year (high luminosity runs). 

Supersymmetry (SUSY) must be broken and the spar- 
ticle mass spectrum depends on the SUSY breaking 
mechanism. Measurement of the sparticle masses pro- 
vides a way to probe the origin of the SUSY breaking in 
nature. The sparticle mass measurement at the LHC has 
been therefore extensively studied. 

Among the sparticles the third generation squarks, 
stops (tj) and sbottoms (pi) (i = 1,2), get special im- 
prints from physics at the very high energy scale. One 
of the examples is found in a model with universal scaler 
masses. Although the scalar masses are universal at a 
high energy scale, the third generation squarks are much 
lighter than the first and second generation squarks due 
to the Yukawa running effect. On the other hand, some 
SUSY breaking models, such as the flavor U(2) model 
[3] or the decoupling solution [4], and the superconfor- 
mal model [5] for the SUSY flavor problem, have non- 
universal boundary conditions for the third generation 
mass parameters at the GUT scale. In addition, the 
tL~t R-Higgs trilinear coupling A t is comparable to the 
gluino mass at the weak scale when the origin of the 
SUSY breaking comes from the GUT scale or the Planck 
scale physics. This appears in the stop mass matrix with 
a large coefficient m t ; m 2 LR = m t (A t — jUcot/3). The 
stop left-right mixing is therefore expected to be sizable, 
leading to an even lighter stop mass compared to those 
of other squarks [6]. It should be stressed that the stop 



masses and the mixing are very important parameters 
to predict the light Higgs mass [7], or the rare B decay 
ratios [8]. 

We may be able to access the nature of the stop and 
sbottom at the LHC provided that they are lighter than 
the gluino (g) . They may copiously arise from the gluino 
decay if the decay is kinematically allowed. The relevant 
decay modes for bi (i = 1,2), t\, to charginos xf (j — 
1,2) or neutralinos Xj U = 1j2, 3, 4) are listed below 
(indices to distinguish a particle and its anti-particle is 
suppressed unless otherwise stated) [21] , 
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In previous literatures [2, 9], the lighter sbottom b\ is 
often studied through the mode (1)2, namely the bbx\ —> 
bbl + l~Xi channel. This mode is important when the 
second lightest neutralino X2 nas substantial branching 
ratios into leptons. The difference of sparticle masses 
such as (rrig — ) is determined by measuring kinemat- 
ical end points of invariant mass distributions for signal 
events. 

In a previous paper [10] we proposed to measure the 
edge position of the m t b distribution for the modes 
(III)i and (IV) 11, where m t b is the invariant mass of 
a top-bottom (tb) system. The decay modes are ex- 
pected to be dominant in the minimal supergravity model 
(MSUGRA), since the branching ratios Br(6i(ii) — > 
t(b)xi) could be as large as 60%. We focused on the re- 
construction of hadronic decays of the top quark, because 
the mtb distribution of the decay makes a clear "edge" 
in this case. The parton level mtb distributions for the 
modes (III), and (IV) are expressed as functions of m§, 
mi , mr , and the chargino mass m~±: dT/dm t b oc mtb, 

and the edge position (end point) of the mtb distribution 
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M t b for the modes (III)j and (IV)»j are written as follows; 




+ \J ( m i ~ ^ m x ± ~ m t) 2 )i m l, - {m~± + m t ) 2 )| (2) 

The measurement of the end points is sensitive to both 
mjj and , provided that m§ and to^± are determined 
from other measurements. The edge height is then closely 
related to cr(pp — > cLY) x Br(III/IV). In some model pa- 
rameters M t b(III)i is very close to M t b(IV)n. When they 
are experimentally indistinguishable, it is convenient to 
define a weighted mean of the end points; 

w BrClIIjMt^III)! + BrClVjnMt^IV)!! 
46 Br(III) + Br(IV)u 

Br(III) ee Br(m)i+Br(III)n+Br(m) 2 i. (3) 

As the final states bbW from the decay chain g — > bh t — > 
6TUii — > bbW\\ (mode (III)»i) could have an irreducible 
contribution to the tb final state, they are included in the 
definition of M™. 

In the previous paper we showed that the weighted 
mean M t ™ was successfully measured from a fit to the ma 
distribution. We also discussed the interpretation of the 
measurement. In this paper we significantly extend our 
study to the branching ratio measurements, identification 
of the stop and sbottom decays into heavier charginos, 
and the interpretation in the MSUGRA model. 

This paper is organized as follows. In Section II, we 
explain our reconstruction and fitting procedure of the 
edge structure in detail. We use a sideband method 
to estimate the background distribution due to misre- 
constructed events, which plays a key role in the edge 
reconstruction. Monte Carlo simulations show that the 
distribution of the signal modes (III) and (IV) after sub- 
tracting the sideband background is very close to the 
parton level distribution. In Section III, we explore the 
MSUGRA parameter space to check if our fitting pro- 
cedure reproduces the end point M™ b and the number of 
events going through the decay modes (III) and (IV). We 
then turn into more delicate issues such as extraction of 
branching ratios. In Section IV, we discuss the stop and 
sbottom decays into the heavier chargino xf by looking 
into the m t b distribution of events with additional lep- 
tons. In Section V, we study the MSUGRA parameter 
region where the decay modes (III) and/or (IV) are open, 
and discuss how to extract the fundamental parameters 
using the m t b distributions. The LHC's potential to ex- 
tract the top polarization arising from the decay g — > tt\ 
is discussed in Section VI. Section VII is devoted to con- 
clusions. 



II. SIMULATION AND RECONSTRUCTION 

In the MSUGRA model the sparticle spectrum is pa- 
rameterized by the universal scalar mass too, the univer- 
sal gaugino mass M , the trilinear coupling of the scalar 
fields A at the GUT scale (Mqut), the ratio of the 
vacuum expectation values tan (5, and sign of the hig- 
gsino mass parameter /i. In order to demonstrate the 
end point reconstruction, we take an MSUGRA point 
with too = 100 GeV, M = 300 GcV, A = -300 GeV, 
tan/3 = 10 and fi > 0. This corresponds to the sample 
point Al in Table I. The masses and mixings of spar- 
ticles are calculated by ISASUSY/ISASUGRA 7.51 [11]. 
Two different event generators, PYTHIA 6.1 [12] and 
HERWIG 6.4 [13], are utilized to generate Monte Carlo 
SUSY events using the masses and mixings. In this sec- 
tion we show results with PYTHIA. The events are then 
passed through a fast detector simulation program for 
the ATLAS experiment, ATLFAST [14]. This program 
performs jet reconstruction in the calorimeters and mo- 
mentum/energy smearing for leptons and photons, giving 
a list of reconstructed jets as well as identified leptons and 
photons. We generate a total of 3 x 10 6 SUSY events, 
which correspond to an integrated luminosity of 120 fb _1 . 
In addition to the SUSY events, ti events are generated 
using PYTHIA as the standard model background. The 
number of generated events amounts to 1.94 x 10 8 , cor- 
responding to an integrated luminosity of 286 fb^ 1 . In 
the plots in this paper, the ti background is not included 
unless otherwise noted. 

In the detector simulation, a lepton is identified if px > 
5 GeV and | rj | < 2.5 for an electron, and pr > 6 GeV and 
1 77 1 < 2.5 for a muon, respectively. A lepton is regarded as 
isolated if it is separated by R > 0.4 from other calorime- 
ter clusters and the transverse calorimeter energy in a 
cone size R = 0.2 around the lept on is less than 1 GeV. 
The cone size is defined as R = \J (A</>) 2 + (Ar/) 2 , where 
4> and r\ are the azimuthal angle and the pseudo-rapidity, 
respectively. 

By default, jets are reconstructed by a cone-based al- 
gorithm with a cone size R = 0.4 in the detector simu- 
lation [22]. After applying the algorithm, jets having a 
transverse energy (Er) more than 10 GeV are kept as 
reconstructed jets. The b and r tagging efficiencies are 
set to be 60% and 50%, respectively. The energy of the 
reconstructed jet is recalibrated according to its jet fla- 
vor using a parameterization optimized to give a proper 
scale of the dijet mass. 

First, we apply the following selection cuts to the sim- 
ulated events for the tb signal: 

1. The missing transverse energy E™ lss is greater than 
200 GeV, where E™ lss is calculated from the re- 
constructed jets, leptons, photons, and unrecon- 
structed calorimeter energies. This calculation is 
performed before the recalibration of jet energies. 

2. The effective mass m e g is greater than 1000 GeV, 
where m e g is the sum of the missing transverse 
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TABLE I: Sparticle masses in GeV and the total SUSY cross 
section (ctsusy) in pb for the parameter points studied in this 
paper. The cross sections are calculated by PYTHIA. 



energy and transverse momentum of reconstructed 
jets, namely m cfi = E™ iss + Y^aiiP'r- 



3. There are two and only two 6-jets with px > 
30 GeV in an event. 

4. Excluding the two 6-jets and those identified as tau- 
jets, the number of remaining reconstructed jets 
with px > 30 GeV and \r}\ < 3.0 should be between 
four and six, inclusive. 

Distributions of the cut variables are shown in Fig. 1. 
The first two cuts are to enhance the SUSY events 
against the standard model background events. The 
other two cuts are to reduce combinatorial background 
(wrong combinations of jets) in the reconstruction of the 
top quark. 

To reconstruct the hadronic decay of the top quark, we 
take the following steps (i)-(iv): 

(i) We first take jet pairs consistent with a hadronic 
W boson decay with a cut on the jet pair invariant 
mass rrijj: \nijj — m\y\ < 15 GeV. The rrijj dis- 
tribution is shown in Fig. 2(a), where the selected 
mass region is marked as W (the W mass region). 
Although fake W pairs dominate the distribution, 
a small bump due to real W bosons can be seen in 
the mass region. 

(ii) The invariant mass of the jet pair and one of the 
6-jets, rribjj, is then calculated. All possible combi- 
nations of jet pairs and 6-jets are tried in an event 
to select the combination which minimizes the dif- 
ference \mt,jj —mt\- The distribution of the selected 
rribjj is shown in Fig. 3(a). 

(iii) The energy and momentum of the jet pair are 
scaled so that rrijj = mw, and the invariant mass 
rribjj is recalculated. The distribution is shown in 
Fig. 3(b). The jet combination is regarded as a top 
quark candidate if \rribjj — m t \ < 30 GeV. 




Number of b-jets 



Number of remaining jets 



FIG. 1: Distributions of (a) missing transverse energy, (b) ef- 
fective mass, (c) number of b-jets, and (d) number of remain- 
ing jets, for 3 x 10 6 SUSY events at the point Al. Accepted 
regions are hatched. 
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FIG. 2: Distributions of the invariant mass of two jets (a) for 
all possible combinations, and (b) for the combination used 
to reconstruct a top candidate. W and A, B indicate the W 
mass region and the W sidebands, respectively. 



As jets are supplied from gluino or squark decays 
and there are several jets in a selected event, events 
with a fake W boson (a jet pair that accidentally has 
rrijj ~ rriw) still dominate the rribjj distribution. The 
contribution of the fake W boson in the W mass region is 
estimated from the events that contain jet pairs with the 
invariant mass in the regions A: \rrijj — {raw — 30 GeV) | < 
15 GeV and B: \rrijj - (m w + 30 GeV)| < 15 GeV. We 
call them "the W sidebands" . The energy and momen- 
tum of the jet pairs are then scaled linearly to be in the 
W mass region \rrijj — mw\ < 15 GeV. The rrijj distri- 
butions before the linear scaling are shown for the W 
region and the W sidebands in Fig. 2(b). The distri- 
bution of the fake top quark candidates is estimated by 
following the same steps (ii)-(iv) for the scaled jet pair 
(see the hatched histogram in Fig. 3(c)). The "true" dis- 
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tribution is obtained by subtracting the background dis- 
tribution estimated by the W sidebands from the original 
distribution (see Fig. 3(d)). The estimation is based on 
an assumption that the background comes from the jets 
without significant correlation with the &-jets. 
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FIG. 4: Distribution of minimum mu for events having a top 
candidate. Open and hatched histograms are for SUSY events 
and tt events, respectively. The loose lepton cut is also shown. 
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FIG. 3: Distributions of rribjj (a) before and (b) after the 
jet energy scaling, (c) background estimated by the W side- 
bands (the hatched histogram) superimposed to the distribu- 
tion (b), and (d) after subtracting the estimated background. 
The mass cut for the top candidate is indicated by arrows. 
Dashed histograms in (c) and (d) show the distributions of 
tt events, where the integrated luminosity of the tt events is 
normalized to that of the SUSY events. 

The tt production is the dominant standard model 
background. The rribjj distributions of the tt events sur- 
viving at this stage are shown in Fig. 3(c) and (d) before 
and after the sideband background subtraction, respec- 
tively. The integrated luminosity of the tt events is nor- 
malized to that of the SUSY events. To reduce the tt 
events, one of the following two lepton cuts may be used, 
depending on the signal/background situation. 

• Loose lepton cut: If there are high-py isolated lep- 
tons, the invariant mass of any high-pT lepton and 
the remaining 6-jet (mu) should be greater than 
150 GeV. The distribution of the minimum mu in 
events with a top candidate is shown in Fig. 4. 

• Hard lepton cut: An event should have no isolated 
leptons. 

The cuts reduce the fraction of the events having a tt 
pair, where one of the top quarks decays leptonically. 
The m^jj distributions with the lepton cuts are shown 
in Fig. 5. The numbers of selected events are listed in 
Table II. Hereafter we use the loose lepton cut. 

For the remaining events, the top candidate is then 
combined with the other 6-jet, which is not used to re- 
construct the top candidate, to calculate the invariant 
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FIG. 5: (a) and (b): Distributions of rribjj with the loose 
lepton cut before and after sideband subtraction, respectively, 
(c) and (d) : Distributions of rribjj with the hard lepton cut 
before and after the sideband subtraction. The meaning of 
the histograms is as same as in Fig. 3(c) and (d). 



mass of the tb system mtb- The distribution is shown 
in Fig. 6(a). However, the expected tb end point is not 
clearly visible in the mtb distribution due to the fake W 
events. Here we can again utilize the W sidebands for 
the background estimation. The estimated background 
distribution is shown in Fig. 6(b), which is obtained by 
averaging distributions from the sidebands A and B. The 
estimated background distribution is subtracted from the 
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59174 


13340 


8168 2764 
(0.138) (0.207) 


Loose lepton cut 


47171 


10487 


5789 1777 
(0.123) (0.169) 


Hard lepton cut 


26915 


6114 


2671 884 
(0.099) (0.145) 



TABLE II: Numbers of events having a top candidate with 
or without the lepton cuts. The integrated luminosity is 
120 fb _1 . The numbers in parentheses are the ratios of the ti 
events to the SUSY events. 



total distribution shown in Fig. 6(c) contains contribu- 
tion from mis-reconstructed events such as other gluino 
decay chains with W and &-jets, or stop/sbottom pair 
productions. 

We fit the total distribution shown in Fig. 6(c), which 
is made with a bin size Am = 10 GeV, by a simple fit- 
ting function described with the end point M^, the edge 
height h per Am bin, and the smearing parameter a orig- 
inated from the jet energy resolution; 
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FIG. 6: (a) Signal mtb distribution for the sample point Al in 
Table I, (b) the estimated background distribution from the 
sideband events, (c) (a) — (b), and (d) the m t b distribution for 
the modes (III) i and (IV) n in Eq. (1), and the decay modes 
(Ill)ii irreducible to the mode (III) i . 



signal distribution in Fig. 6(c). The corrected signal dis- 
tribution (c) shows a better end point structure compared 
to (a). 

By using the information from the event generator, 
one can show the sideband method is indeed reproduces 
the mtb distribution for the signal modes (III) and (IV). 
Fig. 6(d) is the same distribution as Fig. 6(c) but for the 
events which contain either the decay chain (III)i, the de- 
cay chain (IV)n, or the decay chains (IH)ii, which are ir- 
reducible to the chain (III) i ; g — > bbi — > bWt\ — > bbWxf ■ 
Note that if bW has an invariant mass consistent to a 
top quark, the decay is kinematically equivalent to the 
mode (III) i . Fig. 6(d) shows two end points as expected 
(M tb (III)i = 471 GeV and M t6 (IV)n = 420 GeV), 
demonstrating that the sideband method works well. 

Note that the signals from the modes (III) and (IV) in 
Eq. (1) are significant in the total selected events. The 



To reduce the number of free parameters for good conver- 
gence capability of the fit, we set the smearing parameter 
a to be 10 % of the end point Mf£. This assumption is 
based on the dijct mass resolution of the ATLAS detec- 
tor [2] . We assume that the signal distribution is sitting 
on a linearly-decreasing background expressed by a func- 
tion a(m t b — Mfy) + b. The parameters a and b are also 
determined by the fit. 

The fit results, especially the edge height h, depend on 
the mass range used for the fit. We therefore apply an 
iterative fitting procedure to obtain stable results. The 
initial value is determined by fitting a rather wide 
mass range 350 GeV< m t b < 800 GeV (45 histogram 
bins). In the next step, lower 8 bins and higher 25 bins 
with respect to the Mf£ value are used to obtain new fit 
results including the new Mf b l value. The relatively small 
number of lower bins is chosen because the edge height 
is sensitive to the distribution in the mass region near 
the end point. Including too many lower bins is found to 
degrade the sensitivity. The higher bins determines the 
linearly-decreasing background. This step is repeated un- 
til the fit results become stable. The number of iterations 
is typically three or four. If the edge height h has a large 
fit error, the number of lower bins is increased to 9 or 10, 
until the error becomes reasonably small (typically less 
than 15 %). The fit result is shown in Fig. 7. We obtain 
M$ = 455.2 ± 8.2 GeV and h = 271 ± 23 /(10 GeV), 
where the errors are statistical. 



III. THE mtb DISTRIBUTIONS IN THE SUSY 
MODEL POINTS 

In the previous section we discuss the selection of the 
tb events and the extraction of the edge structure in the 
m t b distribution originated from the gluino decays (III)i, 
(Ill)ii or (IV) ji- Two values are obtained from the mtb 
distribution; the end point Mf£ is directly related to the 
stop, sbottom, and gluino masses, while the edge height 
h is sensitive to the gluino production cross section and 
the decay branching ratios. 

This section is aimed to compare the reconstructed val- 
ues with expectations. We will see that the kinematical 
quantities Mf£ and the peak value of the effective mass 
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FIG. 7: A fit to the m t b distribution at the point Al. 
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distributions (M e g) agree very well to M x 
respectively. On the other hand, the reconstructed edge 
height h, which is related to the total number of gluinos 
decaying through the modes (III) or (IV), is dependent 
on the event generator, and also on the decay patterns 
of the neutralino, chargino and squarks. We clarify key 
issues to reduce the uncertainties in this section. 



A. Model parameters 

The gluino decay widths and the decay kinematics de- 
pend on the squark masses and mixings. The charginos 
and neutralinos in the decay chains (I)- (IV) further decay 
into jets/leptons, and the decay patterns depend on the 
electroweak SUSY parameters such as /z, Mi, M2, tan/3 
and the slepton masses. The nitb distribution is there- 
fore a function of all SUSY parameters. To check the 
validity of our reconstruction and subtraction scheme, 
we study the m t b distributions for various sets of SUSY 
parameters. 

The reference point Al introduced in the previous sec- 
tion corresponds to the MSUGRA model with a large 
negative Aq value. The stop and sbottom masses can 
be changed by varying the Ao value without changing 
other sparticle masses drastically (see Section V for the 
detail). We thus make another MSUGRA point A2, by 
changing only the A value from —300 GeV to 300 GeV 
from the point Al. In addition, two non-MSUGRA 
points Tl and T2 are made from the point Al, where the 
chargino/ncutralino sector and the gluino mass are kept 
unchanged but the stop masses and mixing are modified 
by changing the mass parameter m^ R . The stop mass 
parameters and the mixing angle at the low energy scale 
are listed in Table III. 
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TABLE III: The lighter stop mass and the relevant SUSY 
mass parameters in GeV and the stop mixing angles for the 
points Al, A2, Tl and T2. 
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TABLE IV: The universal parameters at the GUT scale for 
the points we study. Units are in GeV except tan (3. 



Next we choose MSUGRA points B, C, G and I, which 
were discussed in the paper [15]. They satisfy nig > irif + 
m t so that the decay mode (III)i is open. The points also 
satisfy the condition m g < 1 TeV, which is required for 
a statistical reason. The gluino production is dominated 
by the process qg — > qg in the case of rriq 
production cross section is parameterized as [16] 



rrig, and the 



<r(q§) + °(?9) = i-74 x ( 



m 5 



uv) 



-5.78 



(pb). (5) 



For example the cross section is 14 pb and 3.1 pb for 
m g — 700 GeV and 900 GeV, respectively. Assuming 
that the branching ratios and the reconstruction effi- 
ciency are equal to those of the point Al, the gluino 
mass must be less than 1 TeV to reconstruct more than 
2000 SUSY tb events after the sideband subtraction for 
an integrated luminosity of 100 fb _1 . 

Finally we select two more points El and E2, where 
only the decay chain (III)i is kinematically open and the 
other two body gluino decays into squarks are closed. 
This happens when niQ and Ao are large. As the SUSY 
events at these points typically contain four bottom 
quarks, the combinatorial background for the tb recon- 
struction is significantly large. 

For the points we study in this paper, we list the rele- 
vant sparticle masses and the universal parameters at the 
GUT scale in Table I and Table IV, respectively. At each 
SUSY point we generate two Monte Carlo data samples, 
each having 3 x 10 6 events, where one is generated by 
PYTHIA and the other by HERWIG. The cross sections 
are summarized in Table I. The SUSY cross section for 
the point C is the smallest among them, where the gener- 
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ated 3x 10 6 events correspond to an integrated luminosity 
of 600 fb" 1 . 



B. Reconstruction of kinematic variables 

The fit results of the mtb distributions are listed in 
Table V, where we follow the fitting procedure described 
in the previous section. The weighted end point M t ™ is 
defined in Eq. (3) and the relevant branching ratios are 
listed in Table VI. In Table V, N a \\ is the number of tb 
events after the sideband subtraction, while N c( i gc is the 
number of tb events after the subtraction, with one and 
only one gluino decaying through the mode (III)i, (IH)ii 
or (IV) ji. Generator information is used to obtain A e d ge . 



The relation between 
Fig. 8. 

weighted end point 



Mi 



The fitted value Mf b l 



tb 



and M™ b is shown in 
increases linearly with the 
in the mass range 350 GcV< 
M^ b < 600 GeV, and tends to be lower than M™. In 
the previous literature [2, 9], the lower mass value was 
understood as the effect of particles missed outside the 
jet cones. In the literature invariant mass distributions 
of the events with same flavor and opposite sign leptons 
are studied, which comes from the cascade decay of the 
squarks q — > qx° — * qll — * Q^Xv The end point of the 
rriju distribution is lower than the parton level qll end 
point by about 10%. 

In the study of the squark cascade decay q — ► qllxl, 
the end point agrees better with the parton level one by 
changing the jet cone size to R — 0.7. In our study at the 
reference point Al, the reconstructed number of events 
is significantly reduced for R = 0.7. We find that using 
the jet cone size R = 0.5 leads to a better reconstruction 
and a larger Mf b , which is closer to M t ™. Such compar- 
ison might be useful to estimate the true end point. We 
discuss the dependence on the jet finding algorithm and 
the cone parameters in Appendix. 

We note that our definition of M™ b might be too sim- 
ple if the reconstruction efficiencies of the decay modes 
(III) and (IV) are very different. In addition, the 
weighted average should not be applied when M t t(III)i 
and M t fc(IV)n differ by more than 80 GeV, as we use the 
events with m tb > Mf b % - 80 GeV for the fit. 

We now discuss the relation between M e g and m g , m q 
using the tb samples. We find that sum of the masses 
m g + niq has a linear dependence on M e g as shown in 
Fig. 9. The deviation of the sample points from a linear 
fit is less than 5%. The plot is for HERWIG samples 
and we find that PYTHIA and HERWIG give consistent 
results for the relation between M e g and the masses. 

The tb sample contains two and only two &-jets origi- 
nating from a gluino decay, therefore the m g dependence 
of M e ff is expected. One may wonder that the gg produc- 
tion might affect the M e g value. This is, however, not 
the case for the points we study. Indeed, the points El 
and E2, where m g » m g , also satisfy the linear relation 
between M e g and rriq + m g . Note that the quark parton 
distribution is harder than the gluon parton distribution. 





gen M^[GeV] M%[GeV] h/{WGeV) 


A'edgc Ngl\ 


Al 


PY 459 455.2 ± 8.3 271.4±22.7 
HW 434.5 ± 5.8 354.8±23.3 


5846 10487 
6685 11470 


A2 


PY 409 442.0 ±17.5 153.0±20.6 
HW 394.4 ±9.5 190.6±21.8 


3064 7525 
3095 7805 


Tl 


PY 468 460.3 ± 5.4 327.0 ± 21.6 
HW 452.0 ± 3.9 447.5 ± 23.5 


6620 11659 
8170 14050 


T2 


PY 429 434.5 ± 8.1 223.2 ± 21.6 
HW 416.6 ± 5.2 321.0 ± 23.2 


4461 8466 
5378 9592 


B 


PY 371 385.6 ± 6.3 226.9 ± 19.6 
HW 361.7 ± 7.3 223.5 ± 21.3 


2801 5396 
3105 5935 


C 


PY 557 548.2 ±14.1 142.4 ± 17.5 
HW 556.3 ±9.4 178.2 ± 17.9 


4026 11228 
4395 12704 


G 


PY 533 498.5 ± 8.6 244.1 ± 23.8 
HW 506.9 ± 6.4 325.5 ± 22.9 


5784 13630 
6248 15039 


I 


PY 507 497.8 ± 7.3 289.7 ± 24.0 
HW 492.9 ± 5.2 383.9 ± 24.5 


6016 13752 
6661 14968 


El 


PY 360 345.5 ± 5.3 270.6 ± 23.7 
HW 348.0 ± 6.3 251.5 ± 23.9 


2778 4595 
3169 5167 


E2 


PY 453 430.8 ± 7.5 352.0 ± 33.9 
HW 444.7 ± 8.0 324.6 ±30.9 


5577 16490 
4394 15724 



TABLE V: Fit results of the mtb distributions and numbers 
of tb events after the sideband subtraction N c dgc and iV a ii for 
3 x 10 6 Monte Carlo events. The fit does not include the ti 
background, and the loose lepton cut is applied. Y a ii is the 
total number of tb events, while Y e d ge is the number of events 
with one and only one gluino decays into mode (HI)i, (III) ji 
or (IV)a. "PY" is for PYTHIA and "HW" is for HERWIG. 



Therefore a(pp — ► qg) cr(pp — > gg) in a wide MSSM 
parameter region. This is why M e s becomes a very good 
function of m g + rriq. 

In Rcf. [17] the relation between M e s and the effective 
SUSY scale is studied, where the effective SUSY scale is 
defined as the cross section weighted mean of the masses 
of two sparticles initially produced in pp collisions. On 
the other hand, we actively select the gq productions by 
requiring two tagged 6-jets. This leads to a clear depen- 
dence of M e ff and rriq ± m g . 



C. Number of tb events 

We now discuss the relation between the edge height 
h, the number of reconstructed tb events, and the recon- 
struction efficiencies. The total number of the "edge" 
events Y d gc arising from the decay chains (III) and (IV) 
may be estimated from Mf b , h per the bin size Am as 
follows, 



N, 



A fit 



h ( m t 
2 \M% 



1 



m t 



Am 



(6) 



This formula is obtained by assuming the parton level 
distribution, and equating the minimum of the mtb dis- 
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(III)i 


(IV)n 


(IV)21 


(HI) n 


(111)21 


sum 


bbX 


Al 


11.0 


6.7 


1.4 


3.4 


2.7 


25.3 


43.4 


A2 


3.1 


6.5 


1.6 


1.4 


0.4 


13.1 


32.0 


Tl 


24.5 


3.2 


0.8 


5.0 


3.0 


36.5 


56.3 


T2 


4.3 


9.9 


2.2 


0.5 


2.1 


19.0 


36.2 


B 


4.1 


8.2 


2.3 


0.9 


1.7 


17.3 


33.5 


C 


7.2 


5.3 


1.3 


0.9 


0.8 


15.4 


38.5 


G 


6.6 


7.5 


1.2 


0.5 


0.8 


16.6 


40.4 


I 


6.2 


11.1 


0.7 


0.0 


0.7 


18.7 


47.3 


El 


78.5 














78.5 


99 


E2 


42.6 














42.6 


98 



TABLE VI: Branching ratios of gluino cascade decays in %. 
The decay modes (III) and (IV) are defined in Eq. (1). The 
"sum" is total of all (III)i, (III);i and (IV) n decay modes. 
The " bbX" is the branching ratio of the gluino decaying into 
t i or bi , so that two bottom quarks appear in the decay prod- 
ucts. 

For 3xl0 6 SUSY events 




300 



300 350 400 450 500 550 600 

M? b [GeV] 

FIG. 8: Relation between M™ b and Mf b for the sample points. 
The solid line corresponds to M™ b = Mf h l and dashed lines to 
M t ™(l ± 0.02) = Mf b \ Bars with a diamond and a circle 
correspond to PYTHIA and HERWIG samples, respectively. 



tribution from the decay chain (III) or (IV) to mt . This is 
a good approximation for reasonable SUSY parameters. 

In Fig. 10(a) and (b) we compare N cdgc in Table V 
with JVfit. We find a very good agreement between them 
both for the PYTHIA and HERWIG samples. On the 
other hand, the correlation between N a \\ and N&t is much 
worse as shown in Fig. 10(c). The number N a \\ receives 
contribution from other gluino cascade decay chains such 
as the modes (I) and (II), as well as contributions from 
the stop and the sbottom pair productions. By the end 
point fit, we extract number of events coming from only 
the modes (III)i, (III) a and (IV) a . 

-^cdgc, and then A^ t , must be related to the number of 
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FIG. 9: Relation between M e g and trig + rriq for the sample 
points in Table I (HERWIG samples) . The line shows a linear 
fit of m g + rriq as a function of M e g. 



gluino decays through the cascade decay chains (III) and 
(IV) via the reconstruction efficiencies. For the points 
we study, the number of produced gg events, N(gg), 
is typically 10% to 14% of the total SUSY production 
events, while the number of q*g and qg production events 
N(q*g) + N(qg) ranges from 42% to 51%. The gluino de- 
cay branchings ratios are listed in Table VI. The number 
of events -/V proc j, where one gluino decays through the 
modes (III) or (IV) and the other squark or gluino decay 
does not produce any bottom quark, is given as follows; 

AWd = 2N(gg)(l-Bv(g^bbX))Bv(edge) 
+ (N(gq)+N(gq*))Bv(edge), 
Br(cdgc) = Br(III)i+Br(III)ii+Br(III) 2 i 

+Br(IV) n , (7) 

where Br(g — > bbX) is the branching ratio of the gluino 
decaying into stop or sbottom, thus having two bottom 
quarks in the final state. The reconstruction efficiency of 
the tb edge mode, etb, is given as 



e tb = N edge /N pmd . 



(8) 



If the efficiency does not strongly depend on uncertainty 
in hadronization and the model parameter dependence 
can be corrected from other measurements, we can ex- 
tract the number N prod from the experimental data. 

The major uncertainty in hadronization may be esti- 
mated by the generator dependence of the reconstruc- 
tion efficiency etb- The edge height h for the HERWIG 
sample is significantly larger than that of the PYTHIA 
sample in Table V, except the points El and E2, and the 
difference is more than 20% at many points. We note 
the difference is small before we apply the sideband sub- 
traction. For example, the numbers of tb events before 
and after the sideband subtraction are 9695 (10180) and 
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FIG. 10: Relations between iV cdgc and jV flt for (a) PYTHIA 
and (b) HERWIG. (c) Relation between N M and jV flt . The 
line in (a) and (b) shows Nm = N a a gc , while the line in (c) is 
a result of linear fit. 



for the HERWIG sample. Indeed, the major difference 
of the two generators are in the fragmentation scheme. 
PYTHIA is based on the string model, while HERWIG 
is based on the QCD model. Currently we are only able 
to say that understanding the nature of fragmentation is 
essential for the interpretation of Na t . We compare the 
HERWIG and PYTHIA samples in Appendix. 

Let us discuss other effects that might change the 
reconstruction efficiency. In Fig. 11 we plot the re- 
construction efficiency e tb as a function of rrig by nor- 
malizing the efficiency to that for the point Al. Here 
^Vprod is calculated using Eq. (7), the generator informa- 
tion on N(gg), N^^g), and the Monte Carlo inputs 
of the gluino branching ratios. The length of the bar 
in the plot shows twice of the statistical error of N c a gc , 
6{e tb /e tb (Al)) = 2(e tb /e tb (Al))/ y/N~^ e . We don't dis- 
cuss the points El and E2 here because events with two 
bottom quarks are not the dominant signature of squark 
and gluino production at these points. 

The efficiencies are very close to one another among 
the points Al, A2, Tl, and Tl. Although the stop mass 
is very different among the points, the mass has little 
influence on the efficiency. Note that these points have 
the same gluino masses, and almost the same parameters 
for the chargino and neutralino sector. For the other 
points, the efficiency can also be expressed as a linear 
function of nig except the sample point C. The gluino 
mass dependence arises from the pre-selection cut m e ff > 
1000 GeV. 

The reconstruction efficiency at the point C (a cir- 
cle in Fig. 11) is lower by 20% from the linear fit (the 
dashed line in Fig. 11). This is due to the loose lepton 
cut which reduces the tt background. As the chargino 
decay branching ratio into leptons is high (about 47%) 
at the point C, the lepton cut kills a significant fraction 
of tb events. The effect is estimated by comparing N e ^ ge 
without any lepton cut and A^dgc with the loose lepton 
cut, at the points C and Al. The efficiency at the point 
C without any lepton cut is plotted in Fig. 11, and is 
close to the linear fit of the other points. 

The gluino mass dependence of the reconstruction effi- 
ciency would be easily corrected by estimating the gluino 
mass from the effective mass distribution. The uncer- 
tainty from the chargino decay branching ratio may be 
corrected by studying tb final states with leptons as well. 
Therefore the parameter dependence of the efficiency 
may be corrected by analyzing the Monte Carlo and real 
data. Efforts to determine the MSSM parameters are of 
course very important to this purpose. 



2462 (2949) for the PYTHIA (HERWIG) sample, respec- 
tively, at the point Al. The number of events before the 
sideband subtraction only differs by 5% between HER- 
WIG and PYTHIA, while the difference increases to be 
more than 20% after the subtraction. This indicates that 
the resolution of the W boson mass is somewhat better 



D. Branching ratios 

In the previous subsection, we find that the reconstruc- 
tion efficiencies significantly differ between the PYTHIA 
and HERWIG samples. However, as we can see in Ta- 
ble V, the generator dependence mostly cancels in the 
ratio Nat/N a \i. Therefore the ratio may play a role to 
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FIG. 11: Reconstruction efficiencies of the model points in 
Table I relative to that of the point Al as a function of m§. 
Point C is plotted twice with/without the loose lepton cut. 



determine the fundamental parameters. For example, 
when the contributions from the stop or sbottom pair 
productions are negligible, Na t /N a \\ is determined by the 
branching ratios to the third generation squarks. How- 
ever, there are several aspects one must consider. 

JV a n consists of contributions from various decay 
modes, even if the contributions from the stop or sbottom 
pair productions are negligible: 

• The ti final state from the decay chains (II) con- 
tributes to JV a n. The branching ratio is typically 
1/3 of that of the tb final state for \i > M 2 > M x . 
The ti final state is reconstructed as tb with a high 
probability, although the rritb distribution does not 
have an edge structure. 

• We expect that events having no t ^ bW de- 
cay would be eliminated by the sideband subtrac- 
tion, and should not contribute to N a \\. However, 
some events actually remain, because of W and Z° 
bosons from the decay of charginos and neutralinos, 
as well as accidental W and Z° bosons from other 
cascade decay chains. An example where x° decays 
dominantly into the Z° boson will be discussed in 
the next subsection. 

In addition, the sideband subtraction is not neces- 
sarily perfect, as will be discussed in Appendix. For 
example, the efficiency of the mode (1)2 is roughly 
a half of that of the mode (III)i after the sideband 
subtraction at the point Al (see Appendix). 

• The cuts to reduce background events may induce 
decay mode dependence of the efficiency. For ex- 



ample the lepton cut could efficiently reduce the 
edge mode (III) if Br(x^ — > IX) is large as we have 
seen in Fig. 11 for point C. 

The interpretation of the ratio N^t/N^w becomes more 
complicated for the following two cases: 

< rrig < rriq, squarks mostly decay into the 



If 



l i(b) 



gluino. As discussed previously, the squark/gluino 
production and decay typically produce events with 
four 6-jets. They are identified as a two &-jet event 
when two &-jets are tagged and the other two are 
mis-tagged. These events decrease the fraction 
from the decay chain (III) and (IV) in the two &-jet 
events. Assuming a tagging efficiency of 60% for a 
single 6-jet, 35% of the four bottom quark events 
are identified as two 6-jet events, and only 1/3 of 
them contain a &-jet pair originated from a single 
gluino decay. The wrong 6-jet pair combinations 
significantly dilute the edge mode. Study of events 
with three &-jets and four &-jets is necessary for this 
case. 

• The tb signal from the stop pair production could 
be significant if the stop is much lighter than the 
gluino. For example, only 4.7% of the SUSY events 
comes from the stop pair production at the point 
Al, while it amounts to be about 17% at the point 
TI. 

In the MSUGRA model, both the right-handed stop 
and the left-handed sbottom is lighter than gluino in a 
broad parameter space. As a result, the decay modes 
which involve W bosons (modes (II), (III) and (IV)) dom- 
inate over the gluino decays to bbX in the region. Because 
the events with W bosons remain after the W sideband 
subtraction, the reconstruction efficiency ttb is expected 
to be similar for those decay modes. Thus, if the con- 
tributions from the stop or sbottom pair productions are 
negligible, the numbers of events with two bottom quarks 
are given approximately as 

iV fit ~ e t6 Br(edge)x 

[2N(gg) (1 - Bi(g -> bbX)) + N(gq) + N(gq*)} , 
AU ~ e tb Br(g ^ bbX)x 

[2N(gg) (1 - Br(~g -> bbX)) + N(gq) + N(gq*)} • 

(9) 

If this simple formula holds, the ratio A^t/^aii 
should provide the ratio of the branching ratios 
Br(edge)/Br(5 -> bbX). 

This is illustrated in Fig. 12, where we plot the ratio 
A^edge/iVaii as a function of Br(cdge)/Br(.g — > bbX). Here 
we plot A^odgc/A^aii instead of Nf^/N^u, because we have 
already seen Nfn <~ N e d ge , and the statistical fluctuation 
of N c( ig C is small (~ 2%) with a help of the generator 
information. Some points in the plots are away from 
the line N cdgc /N a i\ = Br(edge)/Br(g -» bbX): The point 
C is off because the chargino has large branching ratios 
into leptons as discussed earlier. At the point TI, the 
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FIG. 12: Relation between Nedge/N^n and Br(edge)/Br((? — > 
bbX). 

t\t\ productions contributes to -/V a n. The points El and 
E2 are significantly off because the first and the second 
generation squarks dominantly decay into the gluino. 

The systematics due to the chargino and neutralino de- 
cay patterns are difficult to estimate from the LHC data 
only. By combining the data from a future e + e~ linear 
collider (LC) of y/s = 500 ~ 1000 GeV such systematics 
would be significantly reduced. Note that SUSY points 
with m-+ ^ 300 GeV are studied in this paper, which is 
within the reach of a TeV scale LC. 

If events with four 6-jets dominate the SUSY events 
at the LHC, study of the events with four or three b- 
jets would be important, which is out of the scope of this 
paper. If <C mg, it is important to constrain the stop 
mass in order to estimate the contribution from t-\t\, as 
the cross section increases rapidly as decreases. The 
weighted end point M t ™ would provide the information 
on mj . 

E. Snowmass points 

In this subsection we check if we would reconstruct fake 
end points by any chance when the mode of our interest is 
insignificant. To this purpose, we generate a set of Monte 
Carlo data for the Snowmass points [18]. The Snowmass 
points and slopes (SPS) are a set of benchmark points 
and parameter lines in the MSSM parameter space cor- 
responding to different scenarios in the search for SUSY 
at the present and future experiments. Some of those 
points correspond to the case where decay modes (III) or 
(IV) are not dominant, and we should not expect to see 
the edge structure in the m t b distribution. 

We list the fit results at SPS 1-2, and 4-6 in Table VII. 
Here SPS 3 is not listed as it is the same as the point C. 



The M$ is lower by 20 GeV than the M? b at SPS 1 and 
SPS 6, similar to the points discussed in Subsection IIIB 
(see Table V). 

SPS 2 is so-called the "focus point" , where the scalar 
mass at the GUT scale mo is larger than the GUT scale 
gaugino mass Mq. The gluino two body decays into a 
squark and a quark are not kinematically open, therefore 
the dominant decay modes of the gluino are the three 
body decays; Br(g -> tbxf) = 25%, Br(g -> tbxf) = 
20%, and Br(g -> ttxf s) = 20%. The m tb distribu- 
tion should not have any edge structure because it dom- 
inantly consists of the gluino three body decays. The 
reconstructed edge has only 4a in height, and is statisti- 
cally insignificant. The X 2 value of the fit is also rather 
large, \ 2 =44.7 for 30 degrees of freedom [23]. The fitted 
end point M% = 484.9 ± 24.9 GeV may be related to 
the mass difference between the gluino and the heavier 
chargino, m d = 778.6 GeV and m-± = 296.3 GeV. The 
mtb distribution at SPS 2 is compared with that at SPS 1 
in Fig. 13. 

At SPS 5, the stop is so light that the decay t\ — > b\f 
is kinematically closed. The gluino decay branching ra- 
tio into b\ is large (9%), but the b\ subsequently decays 
mostly into TUti(80%). We thus do not find any signifi- 
cant edge. 

SPS 4 is an example that a straightforward applica- 
tion of our reconstruction technique fails to find M t ™. 
The distribution is shown in Fig. 14. The sbottom is 
light because of large tan (3 , and the dominant gluino 
decay mode is g — > bb\ (78%). The sbottom further 
decays into chargino (36%) or the second lightest neu- 
tralino (30%). Decays of the chargino and the second 
lightest neutralino are dominated by X f -> W X \ (100%) 
and %2 ~~ > Z°Xi (98%), respectively. The dominance of 
the decay into Z° or W is common when the decay into 
the Higgs boson or sleptons are kinematically forbidden. 
Some of the decays g — > bb — > bb\2 ar e reconstructed as 
tb events together with additional gauge bosons, result- 
ing in the m t b distribution having a quasi-edge structure. 
The edge is reconstructed by our fit (see Table VII) . Note 
that the end point of the rribb distribution of the decay is 
403 GeV, while the reconstructed m t b distribution of the 
events with the decay (1)2 has a peak around 500 GeV 
and the distribution ends around 580 GeV, close to the 
kinematical limit m~ g — m% <~ 600 GeV. The edge struc- 
ture is identified as that of the mode (I) 2 with a tagged 
Z° — > II. In the mtb distribution, the second edge struc- 
ture from the decay mode (III) or (IV) is weakly seen 
in the distribution. Assuming that two edges exist, we 
find the second end point Mff = 425.4 ± 15.2 GeV and 
h = 349.4 ± 74.4/(10 GeV). 

IV. EVENTS WITH ADDITIONAL TAGS 

At the points given in Table I, the lighter chargino 
and neutralinos (inos) xt > Xi an d X2 are gaugino-like, 
while the heavier inos x h = X21 X3 an d xi are higgsino- 
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FIG. 13: The m tb distributions at (a) SPS 1 and (b) SPS 2. 
The fit curves are also shown. 



SPS 1 
SPS 2 
SPS 4 
SPS 5 
SPS 6 



M? b [GeV] Mf b l [GeV] h/ ( 10 GeV) 



380.8 363.9 ± 4.8 267.3 ± 20.8 

N/A 484.9 ± 24.9 92.2 ± 23.5 

431.0 495.0 ± 6.2 515.7 ± 33.0 

416.2 442.6 ± 42.8 56.6 ± 69.6 

442.3 416.2 ± 6.2 326.4 ± 24.7 



TABLE VII: Fit results for the Snowmass points. Data sam- 
ple at SPS 2 corresponds to 2x 10 6 events, while other samples 
correspond to 3 x 10 6 events. 



like. Unlike the first and second generation squarks, the 
third generation squarks can decay into the heavier inos. 
This is because the third generation squarks couple to 
the higgsinos by the top (or bottom) Yukawa coupling. 
The identification of the cascade decay chain g — > t/b — > 
\ h is a probe for the supersymmetric version of Yukawa 
interactions. In addition, the decay distribution may be 



FIG. 14: The mtb distribution at SPS 4. The fit curve is also 
shown. 



sensitive to the mass difference nig — m^h ~ M 3 — ii. 

In this section we show an example to select such events 
by requiring additional leptons in the final state. We 
study the point C, where the stop and sbottom decay 
into x h with large branching fractions. At this point the 
gluino decay branching ratios are Br(g — > t\) = 15.1% 
and Br(g — > 6i) = 14%. The squarks further decay into 
the higgsino-like inos with branching ratios of Br(6i — ► 
X h ) = 24% and Br(ii -> x h ) = 13%, respectively. The 
stop and sbottom decay branching ratios are comparable 



to those into xf- 

Some of the decay branching ratios of the charginos 
and neutralinos are listed in Table VIII. The higgsino- 
like inos have large branching ratios into the Z° boson. 
The heavier inos also have non-negligible branching ratios 
into triple leptons, because they may also decay into x° 
which further decays into leptons. 

The signature of the heavier inos from stop or sbottom 
is an excess of the Z° boson or three leptons in events 
with two 6-jets. In Fig. 15 we plot the invariant mass of 
the same flavor opposite sign leptons (m«) for the events 
with rif, = 0, 1 and 2, where n& is the number of &-jets 
in an event. In the plots, accidental lepton pair distribu- 
tion estimated with events with different-flavor opposite 
sign leptons are subtracted. All plots show a common 
structure that the invariant mass distribution increases 
toward the edges around 70 GeV and 105 GeV, which 
correspond to the decay chains of x§ ~~ * Ml ~ * MXi an d 
X% — > Ur — > ^Xi, respectively. The Z° peak is also seen 
in the plots. As ri), increases, the Z° peak height becomes 
more significant relative to the edge structure. The ex- 
cess of Z° — > 11 in the m, — 2 sample indicates the Z n 
boson originating primarily from stop or sbottom Quali- 
tatively, the edges in the mu distribution for the rib = 2 
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decay mode 


Br in % 


decay mode 


Br in % 


xt - Wx°2 


29 


x° - WX+ 


59 


xf^z° x t 


26 


xl - z°x ( l 


1.7 


xt - w 


3.0 


X°3 - Z°X°2 


25 




6.8 






xt - Wxt 


59 


X2 - Hi 


9.8 




1.7 


X2 - U~H 


3.0 


x2 - z° x ° 


1.3 


xr - VFx? 


11 


X4 - Jfe 


3.0 


Xf -» lVL 


37 


X4 -» Mfl 


0.9 




8.9 



TABLE VIII: Some of the chargino and neutralino branching 
ratios in % at the point C. 
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FIG. 15: Distributions of m e + e - + m„+„- — m e + M - — m e-n+ 
with n b = (top- left), 1 (top-right), 2 (bottom-left). The 
bottom-right plot shows the distribution for rib = 2 after sub- 
tracting the scaled distribution for n& = by a factor of 0.125. 
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FIG. 16: Same as Fig. 15, but for events with three leptons. 



the third generation squarks to x h is open. 

After establishing the contribution of the decay t/b — > 
X h in the tagged bb3l and bbZ° events, we may recon- 
struct the tb system in the sample. The tb reconstruc- 
tion is similar to that described in Section II, except that 
the distributions are made without the cuts for the effec- 
tive mass and the lepton, because the standard model 
background is expected to be small now. In Fig. 17 the 
mtb distribution with tagged three leptons is shown. The 
number of events after the sideband subtraction is 242.5. 
The mtb distribution is concentrated just below the ex- 
pected end point which is close to m; — m -± ~ 400 GeV, 
although the statistics is low. The distribution of tagged 
Z° — > 11 is similar and the number of events after the 
sideband subtraction is 142. 



sample arc completely subtracted by the distribution for 
the rib = sample scaled by a factor of 0.125, and the 
Z° peak remains (the bottom-right plot of Fig. 15). 

The same analysis can be performed by using events 
with three leptons. In Fig. 16 we plot sum of invariant 
mass distributions; m e + e - — m e + M - — m e - M + for 

all possible combinations of opposite sign lepton pairs. 
In this mode the Z° peak is more significant. The ratio 
of the edge heights in the mu distribution for the n& = 
2 sample to that of the nb = sample is 0.25. _ This 
means that the events originating from the decay b/t —> 
\ h dominate the three lepton events, and the events with 
rib = arise due to misidentified 6-jets. 

Note that these plots at the point C are made for a 
total of 3 x 10 6 SUSY events, which corresponds to an 
integrated luminosity of 600 fb _1 . The statistics for 100 
fb _1 is not significant. However, this procedure should 
be useful for lighter mass spectrum where the decay of 



V. STOP AND SBOTTOM PROPERTIES IN 
THE MSUGRA MODEL 

In this section we interpret our study of the third 
generation squarks at the LHC in a framework of 
the MSUGRA model. The LHC may determine the 
MSUGRA parameter mo, Mo, Aq, and tan/3 through the 
measurement of various distributions. The decay chains 
q -» <7X2 followed by x° -» ^Xi, Z°Xi, 11, or ll X i are espe- 
cially useful if the kinematical end points of the distribu- 
tions are measured. The effective mass measurement of 
the inclusive SUSY signal is also useful to determine the 
absolute SUSY scale. Among these parameters, the A 
parameter is only weakly constrained. This is because 
the left-right mixings of sbottom and stau are almost 
fixed by /utan/3, and A t is insensitive to A as will be 
described below. One of the possible ways to determine 
Aq is to measure the mass difference between squarks 
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FIG. 17: The m t b distributions for events with tagged three 
leptons. Top-left: signal distribution. Top-right: background 
estimated by the sideband events. Bottom- left: distribution 
after the sideband subtractions. 



or sleptons with different flavors, which is caused by the 
rcnormalization group equation (RGE) running between 
the GUT scale and the weak scale. 

The squark-mass matrix is given as follows, 



( ft f~R 



m 
m 



LL 
2 

LL 



rn 



LR 
2 

RR 




where 



= m~ + riif + m| cos 2/3 (- 



rn 



RR 



LR 



• sin Vw ) 



t R ' 



-m| cos 2/3 sin 2 8w, 



[i cot /3) 



(11) 



for stops, and 



m LL = m 2 { l + ml + m 2 z cos 2 [3 (-- + - sin 2 9 w ), 



2 2 2 1 2 no ■ 2 n 

m RR = m l R + m b - cos 2/3 sin 6 W , 
m LR = m b {Ab - At tan /3) 



(12) 



for sbottoms. 

In the MSUGRA, the SUSY-breaking parameters at 
the weak scale in these mass matrices are evaluated by 
the RGEs with the universal GUT scale boundary condi- 
tions. When the bottom- Yukawa coupling is negligible, 
it is convenient to present the SUSY-breaking parame- 
ters by the infrared-fixed point value for the top- Yukawa 
coupling constant J//[19], 



Vf = 



F'(t) 
6F(ty 



(13) 



Here, 
F(t) 



o \acuT ) \aauT J \acuT ) 



where t = l/(47r) log (i 2 /Mq UT and At is the renormaliza- 
tion scale, acur is the gauge coupling constant at the 
GUT scale. The SUSY breaking parameters for the third 
generation at the low energy scale are approximately 
given using the value as 



m,- 

tL 

2 



rn 



In 



0.5ml 
m 2 (l- 



m L 
A t ~ 



0.7M| - 
0.5M? - 

-(0.33A 2 -0.51A Afg)(l-O, 
rog + 0.8Mf, 



- (0.16i4§ - 0.25A M S ){1 - £), 





-0.7M S + A)(1-0. 



A b ~ -1.2M S + 0.8 A 
up to 0((1 -0 2 ). wh ere 



yVvt 2 ) 



203 sin [3 (GeV) 



(15) 



(16) 



Here we use the GUT relation for the gaugino masses, 
Mg ~ 2.8M . For the pole top mass m p t ole = 175 GeV, £, 
ranges 0.7 < ^ < 0.9. 

From Eq. (15), the qualitative behavior of the SUSY- 
breaking parameters is clear. The reduction of the m? 



and 



is more efficient for A^Mg < 0, although the 



dependence is suppressed by an overall factor of 1 — £. 
The right handed stop mass m t - B is also insensitive to 
mo. In a moderate tan/3 region, t\ and b\ is almost 
right-handed and left-handed, respectively. 

In Fig. 18 we show the masses of g, t\, and b\ as func- 
tions of Ao, where the other MSUGRA parameters are 
fixed to be m = 230 GeV, M = 300 GeV, tan/3 = 10, 
and /i > 0. As expected from the discussion above, the 
masses mi and m~ h have sensitivity to Aq. The reduc- 



tion of 



Hi a ^ "% 
I- and m? 

tL t 



by the top- Yukawa coupling con- 



stant is larger when A < 0. 

The stop i\, which is almost Ir, is lighter than the 
gluino unless m is very large. Similarly, since the sbot- 
tom bi is almost 6^, the branching ratio to the edge 
events in Eq. (1), which includes xf > is large since it has 
SU(2) gauge and top- Yukawa interactions. The study of 
the edge events has therefore potentially high feasibility 
in the MSUGRA. Fig. 19 shows the regions, where the 
gluino decays into tl\ and bb\ are open, in a plane of m 
and Aq . In this plot the other parameters are fixed to be 
m g = 707 GeV, tan/3 =5 or 30, and a* > 0. 

Since the top quark mass deviates from the fixed-point 

value, the An dependence on to? and to 2 survives, as 

L R 

shown in Fig. 18. This may lead to a measurable de- 
pendence on Aq of the end point of the m t t distribution. 
In Fig. 20(a), we show M%, M^III^, and M tft (IV) n as 
functions of A . Since b\ is heavy, M t f,(IV) n is lower 
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FIG. 18: Masses for g, t\, and b\ as functions of Ao- Here, 
m = 230 GeV, M = 300 GeV, tan j3 = 30, and /i > 0. 
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FIG. 19: Contours where m g — + m t and m- bi + mt on 

a plane of mo and Ao- The gluino decays into tti and bb\ are 
open below the contours. Here, nig = 707 GeV, tan/3 =5 and 
30, and /i > 0. 



than M t fc(III) 1 . However, Br (g — > tti) is larger than 

Br(g — > bbi) due to the phase space, which makes M t ™ 
closer to M t (,(III) 1 . If is determined with a precision 
of a few % as in our simulation study, the A parameter 
can be evaluated with a precision of ~ 50 GeV. 

In addition to the weighted end point M t ™, the edge 
height h of the m tb distribution is measurable, and the 
height is roughly proportional to the branching ratio of 
the gluino to the edge events arising from the decay 



modes (III) and (IV) in Eq. (1). In Subsection HID we 
proposed to investigate N^t/N a \i where TVgt is the num- 
ber of the edge event estimated by the fit and 7V a ii is the 
total tb events selected. It is related to the normalized 
branching ratio Br(edge)/Br(g — > bbX). The normalized 
branching ratio is sensitive to the branching ratios for i\ 
and 6i, and uncertainty on fragmentation and ^-tagging 
efficiency would be canceled out in the ratio. 

In Fig. 20(b) we show the branching ratio Br(g — > bbX) 
and the branching ratio of the gluino to the edge events 
normalized by Br(g — > bbX). They are decreasing func- 
tions of Ao, since the t\ and b\ is heavier for larger A . 
While the behavior of Br(g — > bbX) is moderate, the 
normalized branching ratio of the edge events is sensi- 
tive to Aq, since it depends on the branching ratios of 
t\ and b\. In Fig. 20(b) the normalized branching ra- 
tio has two kinks, where some decay modes for t\ or b\ 
are open. The first kink around Aq ~ —500 GeV comes 
from a change of the branching ratio of t\. The domi- 
nant decay modes of t\ are b\f and t\i, however, the 
mode t\\ is open around Ao ~ —500 GeV. The decay 
modes t\ — > b\2 an d &i — > t\t are open at the sec- 
ond kink around A ~ —150 GeV. The decay modes are 
not included in the edge events, because the end point 
M t f,(III)2 is significantly lower than M t fc(III)i. The addi- 
tional tags for discussed in the previous section to 
identify the decay works for the region. Since X2 ^ s 
higgsino-like and the Yukawa coupling to b\ is enhanced 
by the top- Yukawa coupling, the decay mode b\ — > txt 
becomes dominant for A > —150 GeV, and then the 
normalized branching ratio become quickly diluted. 

In the MSUGRA, the masses of l\ and b\ are related 
with each other in a broad parameter space. As the deter- 
mination of has been studied in Rcf. [2] , we now dis- 
cuss M™ b and the normalized branching ratio with a fixed 



m; value. We fix mr 

t>i 01 



570 GeV and M § = 707 GeV for 



tan (3 = 5 and 30, and show M™ b (M t t(III)i) as a function 
of mi in Fig. 21(a). Our parameter scan is restricted to 
the range \A \ < 2 TeV. 

We have two disconnected solutions corresponding to 
-1400 (-1850) GeV< A < -280 (470) GeV and 
1350 (1600) GeV < A < 2000 GeV for tan/3 = 5 (30), 
respectively. The A region between the two solutions 
is not allowed because of the experimental constraint of 
f mass or charged LSP. For tan/3=5 (30), m is smaller 
than 460 (750) GeV for A < 0, and 500 (300) GeV for 
A > 0. Since 61 — b L in the MSUGRA, a larger \A Q \ 
corresponds to a larger mo for a fixed , as expected 
from Eq. (15). Nevertheless t\ could be much lighter 
than b\ for a large and negative value of Ao as can be 
seen in Fig. 21(a) [24]. Note that M™ b is determined for a 
fixed t\ mass when to^ is lighter than 370 GeV. If m bi is 
constrained elsewhere, to^ is strongly restricted by the 
M t ™ measurement under the MSUGRA assumption. 

In Fig. 21(b) we show the solution in a M t ™ and 
Br(cdgc)/Br(g — > bbX) plane. The normalized branch- 
ing ratios is almost 1 for M t ™ < 400 GeV (or < 
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FIG. 20: (a) M%, M(III) 1 , and M(IV) n as functions of A . 
(b) Br(g — > bbX) and Br(edge)/Br(<? — > bbX) as functions of 
Aq. Input parameters are the same as in Fig. 18. 
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FIG. 21: For fixed m Sl = 570 GeV and M- g = 707 GeV, (a) 
M? b and m^, and (b^Mfi, (M t6 (III)) and Br(edge)/Br(g -> 
bbX). Here, -2 TeV < A < 2 TeV, tan/3 = 5 and 30, and 
/i > 0. mo is fixed by . The neutralino LSP is assumed. 



VI. MEASUREMENT OF TOP POLARIZATION 



300 GeV). When the mass difference of and m- is 
large, the decay of &i is dominated by Wii. Further- 
more, for m h < 300 GeV, Br(ii -> &Xi ) is 100%. How- 
ever since m= > in this region, squark and gluino 
production goes to the final states having four bottom 
quarks, where we have seen the disagreement between 
the measurement N e d ge /N a n and Br(edge)/Br(g — ► bbX) 
in Section III D. In this case our study must be extended 
to events with more than three tagged 6-jets. 

Note that the decay x° ~~ * I IS mostly closed in Fig. 21. 
The decay is open only at the most right region of 
Fig. 21(b), near the end of the lines of the solutions. 
If this decay is open, the masses of q, I and are model- 
independently determined by the mu, muj and mij dis- 
tributions. When this decay mode is kinematically for- 
bidden, one needs to combine various distributions for 
consistency check of the MSUGRA assumptions. As is 
already shown in Fig. 19, the decay g — ► b/t is open up to 
too ~ m gi providing information of the third generation 
squarks in the wide region of the parameter space. 



Similar to the tau-lepton decay, we may measure the 
polarization of the top quark since it decays to bW via 
the (V — A) interaction. The top quarks from the g, t, and 
b decays are polarized, and the polarization depends on 
the mixing angles for stops, charginos, and neutralinos. 

The bottom quark angular distribution in the polarized 
top quark decay is the following, 

— — ; — —t oc ( — - sin 2 - + 2cos 2 -, (17 
T t dcos6 \m w J 2 2' v ' 

where 6 is the angle between the direction of the bot- 
tom quark and the direction of the top quark spin in the 
rest frame of the top quark. The terms proportional to 
(to(/tow) 2 come from the decay to the longitudinal W 
boson. The bottom quark thus tends to go to the oppo- 
site direction of the top quark spin. 

In the decay mode (III) i , the top quark from a gluino 
decay is polarized to be left-handed (right-handed) if t\ is 
left-handed (right-handed). The polarization is reflected 
on the invariant mass distribution of the bb system (to^), 
if the top quark is relativistic enough in the gluino rest 
frame. Here, one of the bottom quarks comes from the 
top quark decay and the other comes from t\ — > b\ ; . Es- 
pecially, when the invariant mass m t b is close to the end 
point M t b of the decay mode (III) i , the top and bottom 
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FIG. 22: Distribution of m b b in the decay chain (III) i . 
The (dashed) line is for ii = Il{Ir), and 400 GeV< m t b < 
470 GeV. We use the mass spectrum in the sample point Al 
in Table I, and the normalization is arbitrary. 



which go through the decay (III)i . We use the mass spec- 
trum at the reference point Al in Table I. We use events 
with 400 GeV< m tb < 470 GeV to make the distribution, 
and the solid (dotted) line is for the left-handed (right- 
handed) stop. The statistical significance in the differ- 
ence between the left-handed and right-handed stops is 
about 3a for 0(100) events. 

In the above simulation, we neglect contribution of 
other decay chains such as (III)n and (IV) n. The other 
decay chains may contribute to the rribb distribution even 
if we impose that m tb is near the end point. The m bb 
distributions in the modes (III)ii and (IV) n do not de- 
pend on the polarization for the top quark, since the top 
quark distribution in the scalar boson decay is spherical. 
In Fig. 23 we show the mtb distribution around the edge 
region at the reference point Al. The solid line is for 
the total distribution, the dashed line is for the distribu- 
tion excluding the mode (III)i, and the dotted line is for 
the distribution excluding the modes (III)i, (III)n, and 
(III)2i. About a half of the events near the end point 
come from the signal mode (III) 1 , and the ratio between 
the signal events and the rest depends on the MSSM pa- 
rameters. 



VII. CONCLUSIONS 
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FIG. 23: The solid line is for the total mtb distribution at the 
sample point Al, the dashed line is for the mtb distribution 
excluding the mode (IH)i, and the dotted line is for the mtb 
distribution excluding the modes (III) i, (III)n and (III)2i. 



quarks go to the opposite direction to each other in the 
gluino rest frame. Thus, the distribution of the invariant 
mass nibb for events with m t b close to M t f,(III)i is harder 
(softer) for left-handed (right-handed) top quarks. 

In Fig. 22 we show the mbb distribution from the de- 
cay chain (III) i . In this simulation, we use the HERWIG 
generator, since it respects helicities for each particles 
in the processes. We generated a large number of events 



In this paper we study cascade decays g — > 
{ti\ or bbi) — > tb\f at the LHC by reconstructing tb final 
state where the top quark decays hadronically. The m t b 
distribution of the cascade decay has an edge structure. 
The measurement of the end point and the edge height 
of the mtb distribution constrains a combination of the 
masses of g, b, t and % an d the decay branching ratios 
of the particles involved in the decays. 

Through a detailed simulation study, we show in this 
paper the measurement of the end point and edge height 
on a continuum background is indeed possible. Namely, 
the end point of the cascade decay calculated in parton 
level agrees with the reconstructed edge position, and the 
ratio of A^dge (the number of reconstructed edge events) 
and Aan (the number of total reconstructed tb events) is 
understood well by the ratio Br(edge)/Br(g — > bbX). 

The end point and branching ratios depend on the 
mass and left-right mixing of the t and b, as well as 
the chargino and neutralino masses and mixings in the 
MSSM. In the MSUGRA, these sparticle spectrum is ex- 
pressed by a few parameters at the GUT scale. The de- 
cay mode g — > (tti or bbi) — > tbxf is open for a wide 
region of the parameter space where mo < m§. The mtb 
distribution is sensitive to the Aq parameter, the trilin- 
ear coupling at the GUT scale. The distribution is most 
sensitive to A when A ■ M < 0. 

The stop and sbottom could decay both into the heav- 
ier and lighter charginos and neutralinos unlike the first 
and second generation squarks. This is because the third 
generation squarks have the large top (bottom)-Yukawa 
coupling to the higgsinos. A strategy to search for such 
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decays specific for the stop and sbottom is described in 
Section IV. The Yukawa coupling is also related to the 
tL-tR mixing. The polarization of the top quark arising 
from the gluino decay g — > tt depends on the stop left- 
right mixing. The dependence of the rribb distribution on 
the top polarization in the tagged tb sample is discussed 
in Section VI. 

To understand the event distribution better, one needs 
to know the nature of the quark and gluon fragmenta- 
tions into jets. The reconstruction efficiencies are sig- 
nificantly different between the two standard SUSY gen- 
erators HERWIG and PYTHIA, which adopt different 
models for the fragmentation. We point out that the 
smearing of the jet-pair invariant mass originating from 
a W decay affects the reconstruction efficiency. 

The interplay between the LHC and the future LC 
would be useful to reduce the systematics coming from 
the uncertainties of sparticle masses and decay branching 
ratios. They would be reduced dramatically if some of 
the charginos and neutralinos are accessible at the LC. 
Model independent and precise determination of the stop 
and sbottom masses may be possible in such cases and 
will be discussed elsewhere. 
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Appendix: Reliability of Sideband Subtraction 

In this paper, we generate events both by the PYTHIA 
and HERWIG generators. The HERWIG generator uses 
a parton-shower approach for initial- and final- state 
QCD radiations including the color coherence effects and 
the azimuthal correlation both within and between jets 
[13]. The full available phase space for the parton shower 
is restricted to an angular order region, namely, the an- 
gle between the two emitted partons is smaller than that 
of previous branches. On the other hand, the PYTHIA 
generator adopts the string model [12]. The two genera- 
tors predict different tb reconstruction efficiencies, which 
may be considered as the uncertainty in fragmentation. 

Definition of jets also affects the number of recon- 
structed tb events and the reconstruction efficiency. We 
try two algorithms, a cone-based algorithm and a Kt al- 
gorithm (Montreal version), available in the JET Finder 



gen jet 


Mf^GeV] ft/(10GeV) 


PY cone 0.4 
0.5 
0.6 
K T 0.4 
0.5 
0.6 


455.2 ± 8.2 271.4 ± 22.7 
436.1 ± 7.7 272.9 ± 33.4 

461.1 ± 10.4 217.8 ± 22.8 

442.0 ± 4.7 321.3 ± 22.3 

452.2 ± 6.0 305.3 ± 21.3 

459.1 ± 6.1 241.6 ± 31.6 


HW cone 0.4 
0.5 
0.6 
K T 0.4 
0.5 
0.6 


434.5 ± 5.8 354.8 ± 23.3 
460.2 ± 4.9 349.2 ± 22.8 
440.9 ± 7.1 305.3 ± 33.7 
434.9 ± 4.3 406.5 ± 22.1 
460.0 ± 5.5 379.6 ± 33.5 
468.4 ± 5.8 314.3 ± 20.7 



TABLE IX: Fit results at the point Al. We use the cone- 
based algorithm or the Kt algorithm with cone sizes of R = 
0.4, 0.5, 0.6. The fit is based on 3 x 10 6 events. The M? h is 
459 GeV. 

Library [20], which is interfaced to the ATLFAST pack- 
ages. 

The cone-based algorithm merges the calorimeter cells 
around a high Et cell in a fixed cone size R. On the 
other hand, in the Kt algorithm, a cell called a protojet 
i grows until there are no more protojets j which satisfy 
(Ar/) 2 + (A(j>) 2 < (E Ti /mm(E Ti ,E Tj ))R 2 , therefore the 
cone size depends on Ex's, and the shape of the jet de- 
pends on the distribution of the particles in the jet. The 
K T algorithm is based on QCD, and it is considered to 
be advantageous to merge soft jets from an initial parton, 
although experimentally challenging. 

The fit results with different generators and jet finding 
algorithms at the point Al is summarized in Table IX. 
The end point Mf£ becomes closer to M™ b as R increases. 
This is expected because more soft jets are merged to 
leading jets as R increases. The edge height h is larger 
for HERWIG with the K T algorithm. The difference of 
the height for R < 0.5 is more than 25% in Table IX. 
The height significantly decreases for R = 0.6. This is 
because we have to reconstruct a jet pair with relatively 
small invariant mass <~ 80 GeV. 

To see the origin of the difference more closely, we 
compare the rrijj distributions consistent to the top in- 
terpretation in Fig. 24, where, the solid histograms are 
the invariant mass distributions for jet pairs which satisfy 
\m,jj—mw\ < 15 GeV and \rribjj — m t \ < 30 GeV. Dashed 
histograms show the distribution of accidental jet pairs 
in the W mass region, that is, the "fake IF" background, 
estimated by the sideband method. The left plots are the 
distributions of the events arising from the decay chain 
(III) i , while the right plots are the m t b distribution of all 
selected tb events. The distribution at the bottom-left 
(HERWIG and the Kt algorithm for R = 0.5 ) are more 
concentrated around mjj ~ 80 GeV compared to the oth- 
ers (PYTHIA and the cone-based or the Kt algorithm), 
corresponding to better reconstruction efficiencies. 
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TABLE X: Numbers of events before and after the sideband 
subtraction for the events with a gluino that decays through 



mode (III) 1 g — > tti — > tb\f and mode (1)2 g — ► bbi 



bbxl 



sideband subtraction, respectively. Comparing these ra- 
tios to that of the branching ratios Br(I) 2 /Br(III) 1 =0.32, 
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Incompleteness of the sideband subtraction is seen in 
the same figure. The sideband subtraction uses the jet 
pairs with invariant mass 35 GeV< rrijj < 65 GeV and 
95 GeV< rrijj < 125 GeV to estimate the "fake W" back- 
ground. On the other hand, as the jets are defined to 
have >10 GeV, a jet pair near rrijj <~ 35 GeV is sup- 
pressed. We see an underestimate of background events 
in the range 65 GeV< rrijj < 70 GeV, where the dotted 
line and solid line differ significantly. We find that events 
with rrijj < 70 GeV do not contribute to the edge struc- 
ture of the m t b distribution. The events in this region 
might have to be removed from the candidate events. 

The incompleteness of the sideband subtraction can 
also be seen by comparing the reconstruction efficiency 
of the events from specific decay chains. Here we take 
two decay chains for comparison; 
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The decay chain (III)i is relevant to the tb events. On 
the other hand, as the decay chain (I) 2 has no W, the 
sideband subtraction should remove the events arising 
from the decay chain completely, unless there are acci- 
dental weak bosons, or the understimation of the misre- 
constructed events that has just been discussed above. 

The branching ratio Br(x 2 ^ Z°Xi) 1S small (1.5%) at 
the point Al. The decay branching ratios at the point 
Al are 



Br(III)i = Br(.t) -> tii)Br(fi - 

= 0.078 

Br(I) 2 = Br(£ -» 66i)Br(6i 

= 0.025, 
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FIG. 24: Distributions of rrijj consistent to the W boson 
mass under the cuts described in Section II. The MC data is 
generated for the point Al. In the left, we plot the events 
that contain the decay chain g — ► tti — > tbxt only, while 
we plots all events in the right. Dotted lines show estimated 
background. Top figures: PYTHIA with the cone-based al- 
gorithm (R = 0.4). Middle: PYTHIA with the K T algo- 
rithm (R = 0.5). Bottom: HERWIG with the K T algorithm 
(R = 0.5). 
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where the factor 0.7 is the hadronic branching ratio of 
the W boson. 

The numbers of reconstructed tb events before and af- 
ter the sideband subtraction are listed in Table X. For 
the Kt algorithm with R = 0.5, the ratio of the num- 
bers N(I) 2 /N(III)i are 0.229 and 0.15 before and after the 



we see that the requirement \rribjj — m t \ < 30 GeV re- 
duces the contribution of (I) 2 to about 70% relative to 
(III) i and it is further reduced to less than 50% after 
the sideband subtraction. This ratio may be improved 
further by cutting the events with rrijj < 70 GeV or by 
using more sophisticated background estimation. 
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